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1 Abstract 



The SELEX experiment (E781) at Fermilab has observed two statisticaUy compeUing high 
mass states near 3.6 GeV/c^, decaying to A^K^vr^ and A+K~7r"'"7r"'". These final states are 
Cabibbo-aUowed decay modes of doubly charmed baryons S+ and respectively. The 

masses are in the range expected from theoretical considerations, but the spectroscopy is 
surprising. SELEX also has weaker preliminary evidence for a state near 3.8 GeV/c^, a high 
mass state decaying to A+K~7r"'"7r"'", possibly an excited 5++ (ecu*). Data are presented and 
discussed. 



2 Introduction 

The quantum chromodynamics hadron spectrum includes doubly charmed baryons 
(DCBs): S+ (ccd), S++(ccm), and f2+ (ccs), as well as the triply charmed ^2++ (ccc). A 1996 
DCB review [1] collected information on masses, lifetimes, internal structure, production cross 
sections, decay modes, branching ratios, yields, and experimental requirements for optimizing 
the signal and minimizing the backgrounds. DCB works published since then are given in 
Refs. [2, 3, 4, 5, 6, 7, 8, 9, 10]. The doubly and triply charmed baryons provide a new window 
for understanding the structure of all baryons. As pointed out by Bjorken [11], one should 
strive to study the triply charmed (ccc) baryon. Its excitation spectrum, including several 
narrow levels above the ground state, should be closer to the perturbative regime than is the 
case for the DCBs. The (ccq) studies are a valuable prelude to such (ccc) efforts. 

Hadron structures with size scales much less than l/A^cd shouki be well described by 
perturbative QCD. The tightly bound color antitriplet (cc)3 diquark in (ccq) may satisfy this 
condition. But the DCB radius may be large, if it is dominated by the low mass q orbiting 
the tightly bound (cc) pair. The study of such configurations and their weak decays can help 
to set constraints on models of quark-quark forces [12, 13]. Stong [14] emphasized how the 
QQq excitation spectra can be used to phenomenologically determine the QQ potential, to 
complement the approach taken for QQ quarkonium interactions. 

Savage and Wise [15] discussed the (ccq) excitation spectrum for the q degree of freedom 
(with the (cc) in its ground state) via the analogy to the spectrum of Qq mesons, where the 
(cc) pair plays the role of the heavy Q antiquark. Fleck and Richard [12] calculated excitation 
spectra and other properties of (ccq) baryons for a variety of potential and bag models, 
which describe successfully known hadrons. In contrast to heavy mesons, the descriptions of 
light quark (qqq) and singly charmed (cqq) baryons are less successful. We need to better 
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understand how the proton and other baryons are buih from quarks. The investigation of 
the (ccq) system should help put constraints on baryon models, including light quark (qqq) 
and singly charmed (cqq) baryons, since the (ccq) has a quark structure intermediate between 
(qqq) proton and Qq meson structures. 

In the double-charm system, there have been many predictions for the masses of the J=l/2 
states and the J=3/2 hyperfine excitations [9]. Most results are consistent with expectations 
of a ground state mean mass around 3.6 GeV/c^. The (cc) color antitriplet diquark has spin 
S=l. The spin of the third quark is cither parallel (J=3/2) or anti-parallel (J=l/2) to the 
diquark. For (ccq), the J=l/2 states are expected to be lower than the J=3/2 states by around 
80 MeV/c^ [9, 10, 12, 16]. 

Bjorken [11] and also Fleck and Richard [12] suggest that internal W exchange diagrams 
in the S+ decay could reduce its lifetime to around 100/s, roughly half the lifetime of the 
A^. Considering possible constructive interference between the W-exchange and two c-quark 
decay amplitudes, it is possible that this state should have an even shorter lifetime. 

We describe qualitatively the pcrturbativc production mechanism for DCBs. One must 
produce two c quarks (and associated antiquarks), and they must join to a tightly bound, 
small size anti-triplet pair. The pair then joins a light quark to produce the final (ccq). The 
two c-quarks may be produced (initial state) with a range of separations and relative momenta 
(up to say tens of GeV/c). In the final state, if they are tightly bound in a small size (cc) 
pair, they should have relative momentum lower than roughly 1 GeV/c. The overlap integral 
between initial and final states determines the probability for the (cc)-q fusion process. Such 
cross sections may be smaller by as much as 10^^-10^^ compared to single-charm production. 
Aoki et al. [17] reported a low statistics measurement at ^J1 =26 GeVjc^ for the ratio of 
double to single open charm pair production, of 10~^. This DDDD to DD cross section ratio 
was for all central and diffractive events. This high ratio is encouraging for (ccq) searches. 
Cross section guestimates are given in Ref. [1]. 

Brodsky and Vogt [18] suggested that there may be significant intrinsic charm (IC) cc com- 
ponents in hadron wave functions, and therefore also cccc components. The double intrinsic 
charm component can lead to (ccq) production, as the (cc) pairs pre-exist in the incident 
hadron. Intrinsic charm (ccq) production, with its expected high X/ distribution, would 
therefore be especially attractive. When a double charm IC state is freed in a soft coUision, 
the charm quarks should also have approximately the same velocity as the valence quark. 
Thus, coalescence into a (ccq) state is likely. Cross section guestimates are given in Ref. [1]. 

The semileptonic and nonleptonic branching ratios of (ccq) baryons were estimated by 
Bjorken [11] in 1986. He uses a statistical approach to assign probabilities to different decay 
modes. He first considers the most significant particles in a decay, those that carry baryon or 
strangeness number. Pions are then added according to a Poisson distribution. The Bjorken 
method and other approaches for charm baryon decay modes are described by Klein [19]. For 
the S^"*", Bjorken [11] estimated the A'^tt'^K'tt'^ final state to have 5% branching ratio; while 
for the S+ , he estimated the A'^n^K" final state to have 3% branching ratio. One expects [1] 
that roughly 80% of the (ccq) decays are hadronic, with as many as one-third of these leading 
to final states with all charged hadrons. 

3 Observation of DCBs at Fermilab SELEX 

The SELEX experiment (E781) at Fermilab [20, 21, 22, 23] carries out DCB data analysis 
based on the sample of 1630 A+ and 10210 D-meson events from hfetime studies of these 
particles. SELEX used 600 GcV/c beams of 7r~, E~, and protons to produce charm particles 
in Carbon and Copper targets, detecting them with a magnetic spectrometer with high mass 
resolution (pa 5 MeV/c^) and high proper lifetime resolution (pa 20 fs.), and Cerenkov particle 
identification of particles above 25 GeV/c. 
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Figure 1: Schematic of S+ —>■ K~7r"'"A+ 



SELEX has two statistically-compeUing new high- mass states near 3.5 GeV/c^ decaying to 
A+K-7r+ (3460 MeV/c^) and A+K-7r+7r+ (3520 McV/c^) decay channels of DCBs S+ (ccd) 
and (ecu) [24, 25]. They appear to be members of the DCB family of new particles: 
(ccd) DCBs composed of two charm quarks and one down quark; and (ecu) DCBs composed 
of two charm quarks and one up quark. SELEX also has evidence [26] for a (ecu) baryon 
near 3780 MeV/c^, a high mass state decaying to A+K~7r"'"7r"'", which may be a decay mode 
of an excited DCB 5++ (ecu*). The 3520 MeV/c^ state satisfies all expectations for being 
a H+ state. Its mass is compatible with most model calculations for the DCB ground state. 
Its hfetime is shorter than 30 fs at 90% confidence. The 3460 MeV/c^ state has the decay 
characteristics of a state, and a comparable lifetime. However, it is difficult to understand 
the 60 MeV/c^ mass difference between the Q=l and Q=2 states if they are members of the 
ground state isodoublet. It is not yet confirmed that J = 1/2 for both states. 

SELEX analysis and simulations continue for the charmed Lambda (A^) decay-mode DCB 
data. SELEX also searches for the states observed to date in decay channels where the final 
state has a charmed D-meson rather than a Ac-baryon. Continuing SELEX data analysis 
aims to achieve a more complete understanding of the double charm sector, production with 
different beam particles, decay modes, cross sections and their A-dependence, Xp and pt 
distributions. 



3.1 Features of the Selex spectrometer 

SELEX is a 3-stage magnetic spectrometer [20]. The negative 600 GeV/c Fermilab Hyperon 
Beam had about equal fluxes of 7r~ and E~. The positive beam was 92% protons. For charm 
momenta in a range of 100-500 Gev/c, mass resolution is constant and primary (secondary) 
vertex resolution is typically 270 (560) fim. A RICH detector labelled all particles above 
25 GeV/c, greatly reducing background in charm analyses. [21]. The details of single-charm 
analyses involving Aj" — > pK'n'^ reconstructions can be found in [22, 23]. 

3.2 Double-charm Analysis for Ac Decay Modes 

The double-charm search began with the sample of 1630 A+ events (cud) used in the lifetime 
analysis. We ask for a weak-decay vertex lying between the primary vertex and the observed 
A+ decay vertex. A Cabibbo-allowed ccd decay can give a final-state A+, a K~, and a tt"*". 
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Figure 2: (a) The A+fr^7r+ mass distribution in 5 MeV/c^ bins. The shaded region 3.400- 
3.640 GeV/c^ contains the signal peak and is shown in more detail in (c). (b) The wrong-sign 
combination A'^K'^n' mass distribution in 5 MeV/c^ bins, (c) The signal (shaded) region 
(22 events) and sideband mass regions (140 events) in 2.5 MeV/c^ bins. The fit is a Gaussian 
plus linear background. 



The c decays weakly, for example by ccd csu + uu + dd; corresponding for example to a 
ccd A+K~7r+ final state. The event topology contains two secondary vertices. In the first, a 
A;t and a K~7t^ meson pair are produced. This vertex may be distinguished from the primary 
vertex, if the (ccd) lifetime is sufficiently long. The now propagates some distance, and 
decays at the next vertex. The experiment must identify the two secondary vertices. The idea 
is shown schematically in Fig. 1. 

We reconstruct total electric charge states Q=l and Q=2 in separate reconstructions. The 
Q=l sample consists of two oppositely-charged tracks along with a A^. The Q=2 sample 
consists of a positively-charged triplet along with the A+. Because the meson-pair tracks from 
the intermediate vertex may have energies below the RICH detection threshold (25 GeV/c 
for pions), or since they miss the RICH when they are emitted at wide angles, few achieve 
particle identification in the RICH. To be consistent with a ccd — > A+K~7r+ final state, the 
negative track should be a Kaon. We build "right-sign" and "wrong-sign" samples, based on 
calling the negative track a Kaon or a pion, respectively. 

Event selection cuts used here were taken without change from previous single-charm 
studies. For short-lived states, Li/cxi > 1 and the A+ momentum vector must point back 
to the primary vertex within a cut. We have varied the cuts and observe that no signal 
significance depends critically on any cut value. The signals seen here are stable. 



3.2.1 Q = 1 reconstruction, 3520 MeV/c^ 

Fig. 2. shows the K^7r+A+ mass distribution for single-charged baryons. Fig. 2(a) shows a 
S+ candidate at 3520 MeV/c^, consistent with most model calculations. The peak is 3 ± 1 
MeV/c^ wide, narrower than but consistent with simulation (5 ± .5 MeV/c^). The final state 
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Figure 3: — > K~7r+7r+A+ mass distribution in 7.5 MeV/c^ bins. 

meets the double-charm Cabibbo-favored criteria: final state single charm, final state baryon, 
final state K~. The general agreement between right-sign (negative track is Kaon) and wrong- 
sign (negative track is pion) average levels and fiuctuations in Fig. 2 confirms that most events 
are combinatoric background. The signal channel shows a 22-event excess over a background 
of 6.1 ± 0.51 events, for a single-bin significance of 6.3 a. Treating this with gaussian statistics, 
the probability of such an excess is less than 10~^ for a single bin. We searched for a peak in 
the interval 3.2-4.3 GeV/c^, or 110 bins. In each bin there is a Poisson-distributed number 
of events along with a Gaussian fiuctuation on the background. Summing these effects gives 
the overall probability that our search would find such a fiuctuation anywhere in the search 
interval of less than 10~^. 

We calculated the lifetime likelihood for signal and sideband regions in reduced proper 
time t* = M{L — Lmin)/pc for several different lifetime values. In this analysis Lmin = cxi, the 
error on the vertex separation Li in the schematic from Fig. 1. From SELEX single-charm 
lifetime studies, we know the lifetime resolution is about 20 fs. Here, the sidebands show a 
lifetime of 30li\l fs. The signal region result is shorter, with lifetime less than 30 fs at 90% 
confidence. These two results are, of course, statistically consistent with a common short 
lifetime for signal and sideband regions. 

3.2.2 Q=2 Reconstruction, 3460 MeV/c^ 

Fig. 3 shows a candidate at 3460 MeV/c^. The peak is narrow, 6 ± 1 MeV/c^, matching 
our simulation width. The final state meets the double-charm Cabibbo-favored criteria: final 
state single charm, final state baryon, final state K~. The Q=2 sample distribution is quite 
different from that for Q=l. The kinematic threshold for this channel opens near 3.0 GeV/c^, 
yet we have only 2 events below the 3460 mass peak. The background under the peak is clearly 
Poisson-distributed. Making a likelihood fit to a linear background function, we find 9 events 
in the peak, compared to an expected background of 1.6 events. The Poisson probability that 
there is an excess of 7.4 events or more is 5 x 10~^. This search looked for an isospin partner 
to the S+(3520) and was limited to the range 3.3-3.6 GeV/c^, containing 13 bins of width 
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Figure 4: S^"*" — >■ K 7r+7r+A+ mass distribution in 7.5 MeV/c^ bins. The linear background 
is determined from a liklihood fit in the mass range 3.41-3.99 GeV/c^. 

22.5 MeV/c^. The probabihty for a statistical fluctuation up to a 7.4-event excess in any of 
these 13 bins is less than 7 x 10"''. 

There are too few events to attempt a lifetime analysis, either for signal or background. 
We have compared the rest-frame angular decay variables to a phase-space simulation. The 
comparison is good for the signal region. The background tends to be more peaked along the 
flight direction of the candidate, as one might expect for accidental reconstructions. 

One can ask why there are so few candidates compared to S+(3520) candidates. 
Simulation shows that the acceptance is substantially lower for the Q=2 system at masses 
below 3500 MeV/c^ than for a corresponding Q=l system. The SELEX spectrometer cuts off 
the slow pion acceptance from the higher-multiplicity state. This can be seen in the data by 
the presence of low-mass background in Fig. 2 and its absence in Fig. 3. 

3.2.3 Q=2 Reconstruction, 3780 MeV/c^ 

We also searched at higher mass for A+K^7r+7r+ decay of (ecu) [24]. We have evidence 
[26] for a ecu baryon near 3780 MeV/c^, a state which may be an excited DCB S^"*" 
(ecu*). The data are shown in Fig. 4, with a linear fit to the background. Other fits are 
possible, with higher background levels, as in Fig. 5. The background curves are determined 
from the spectrum itself, from background simulation studies, and taking into account also 
the wrong-sign mass plot. The higher background in Fig. 5 is based on a simulation study, 
which shows that the excess near 3660 MeV may be due to combinatorial background of a 
real ccd with a 7r+ from the primary vertex. This leads us to fit from 3650 MeV, to account 
for this effect. The number of signal events drops, the background increases, and the signal 
significance drops from 5. 60" to A.Oa. The width of the peak is roughly three times wider than 
the low lying DCB peaks, which would be possible for a state which decays strongly. For this 
ecu* however, the data do not have a plausible explanation, in that roughly 50% of the signal 
events above background decay weakly (to A+K^7r+7r+) and 50% decay strongly (to vr+S^). 
This is not possible for a single state. We can not now draw unambiguous conclusions. We 
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Figure 5: H^"*" — >■ K 7r+7r+A+ mass distribution in 7.5 MeV/c^ bins. The linear background 
is determined from a liklihood fit in the mass range 3.65-3.99 GeV/c^. 



plan further study of these data and backgrounds. 



3.3 Production 

Production characteristics of the 22 signal plus background events at 3520 MeV/c^ are indis- 
tinguishable from the single-charm sample [23]. The mean pt is 1 GeV/c and mean xp ~ 
0.33. This value of a; is significantly higher than the value xp ~ 0.10 that one expects [1] if 
the (ccd) is produced near threshold in a central collision. 

We compared production of the and states by different beam hadrons. The 
DCBs we see are produced solely by the baryon beams in SELEX data. There are no signal 
candidates from the pion beam. We also checked if there is a dependence on the target nucleus. 
We found [24] that the diamond/copper ratio of the signal events is similar to the sideband 
events, which in turn behave like single-charm production. 

Count rate estimates for DCBs at were given previously [1] based on ccq cross section of 
at most 25 nb per nucleon (25 nb/N). This value is based on a reduction factor of 1000 per 
produced charm quark, considering that the open charm production cross section (25/ib/N) 
is 1000 times lower than the total inelastic cross section 25 mb/N. But SELEX DCB cross 
sections are several orders of magnitude higher, as approximately 60 DCBs are observed in a 
sample of 1630 Ac. Simulation studies suggest that the double-charm states may account for 
as much as 40% of the A^ sample seen in this experiment, a surprisingly high fraction. Work 
is in progress to calibrate the SELEX trigger efficiency factor, in order to fix the absolute 
value of the DCB cross sections. The high DCB cross section is reminiscent of the discovery 
of the S+ baryon in the WA62 experiment at CERN, using a 135 GeV/c hyperon beam [27]. 
That cross section was also anomalously large and still has no good theoretical explanation. 
The FOCUS photoproduction experiment at Fermilab has looked for these states using their 
A+ events. They see very few candidate events and no signal peaks [28]. 
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Figure 6: — > K 7r+7r+A+ mass distribution in 7.5 MeV/c^ bins, with angle cut. Signal 
events are shaded. Wrong-sign background is shown as open histogram boxes. 

3.4 Background Reduction Studies for Ac Decay Channels 

The cleanest SELEX DCB data are for ccd and ecu lowest peaks, while for ecu* the back- 
grounds are significant. Consider the ecu at 3460 MeV. Events outside the signal region show 
a strong preference for the center-of-mass (CM) angle of the negative track to be near 180 
degrees, where the angle is with respect to the A+ direction. This suggests that these tracks 
represent accidental overlaps of a low energy primary track and a real secondary vertex. In 
the CM system, a low energy primary track is interpreted as a back-angle track from the 
secondary vertex. 

Simulation indicates that a cut to remove such events should have little effect on the signal 
region for a phase-space distribution. That is indeed the case in the data. For the ecu at 3460 
MeV, with a cut cos(^^) > —.6, we obtain 1 background and 8 signal events. The spectrum is 
shown in Fig. 6. This angle cut gives result cleaner than without the angle cut, where we had 
previously 1.6 background and 7.4 signal events. Without this cut, the Poisson probability 
that there is an excess of 7.4 events in any of the 13 bins between 3.3-3.6 GeV/c^ is less than 
7 X 10~"^, and the signal significance is 5.8 a. With this cut, the Poisson probability that there 
is an excess of 8 events in any of these 13 bins is reduced to 10~^, and the signal significance is 
7.9 a. We will make similar studies for the ecu* data. The CM angle cut does affect the ecu* 
region strongly, but not in a way that is simple to understand. For ecu, ecu*, SELEX studies 
to what extent the decay angular distribution of the K~ in the rest frame of the K~7r"'"Aj 
system is compatible with a phase-space simulation, both for signal and for background. With 
further simulations, SELEX will try to better understand such experimental distributions, and 
will explore other methods of reducing the backgrounds. 

3.5 Search for double-charm baryons in D-meson decay modes 

SELEX started complementary DCB data analysis of these same states based on our indepen- 
dent data for D-meson decays. Some possible decay modes of interest are: {D~^K~p), {D^K^p) 
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for (ccd); (D+A7r+), (L'0A7r+7r+), {D+K^p), (L'*+A7r+), {D+K-p7r+), for (ecu). This work is 
in progress. 

4 Planning for the CERN COMPASS experiment 

The SELEX efforts will help in the planning of a future state-of-the-art double charm mea- 
surement via the CERN COMPASS experiment, and elsewhere. The Sept. 2002 COMPASS- 
Future meeting at CERN included already a presentation of SELEX data, and initial COM- 
PASS plans in response to consultations with SELEX on the SELEX DCB data [30] . Based on 
the SELEX results, one could expect from COMPASS [1, 30] up to 10,000 reconstructed DCBs 
for a 100 day run. The SELEX data have stimulated more detailed plans/prospects/efforts 
at COMPASS to study DCB and even triply charmed baryon spectroscopy. 

5 Summary 

We reviewed DCB theory. We presented SELEX data for two statistically-compelling new 
high-mass states that decay into a final state A+, K~ and one or two tt"*", as expected for 
double-charm baryon decays. The 3520 MeV/c^ state satisfies all expectations for being a 
state. Its mass is compatible with most model calculations for the double-charm baryon 
ground state. Its lifetime is shorter than 30 fs at 90% confidence. The 3460 MeV/c^ state 
has the decay characteristics of a 2++ state. However, it is difficult to understand the 60 
MeV/c^ mass difference between the Q=l and Q=2 states if they are members of the ground 
state isodoublet. We have not yet tried to confirm that J = 1/2 for both states. We also 
showed evidence [26] for a (ecu) baryon near 3780 MeV/c^, a high mass state decaying 
to A+K^TT+TT"*", which may be a decay mode of an excited DCB S++ (ecu*). The SELEX data 
have been used already in the planning of the CERN COMPASS DCB experiment. Internet 
hnks to SELEX and COMPASS DCB presentations are given in Ref. [31]. 
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